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THE POLAR UPPER ATMOSPHERE AND AURORA 


Kaichi Maeda 

Goddard Space Flight Center, NASA 
Greenbelt, Maryland 

George Ishikawa 

Meteorological Research Institute 
Suginami, Tokyo, Japan 


ABSTRACT 

This paper reviews the following three topics in the physics of 
the polar upper atmosphere and their latest developments: 
(i) Winter anomalies of the polar ionosphere and their interpreta- 
tions, (ii) Aurorae and some related phenomena, and (iii) Noctilucent 
clouds. In each topic, there are still unsolved problems and con- 
troversies that could be settled by direct observations by means of 
rockets and artificial satellites. The importance of rocket-soundings 
in the altitude range between two limits, i.e., the one around 40 km 
for balloons and the other around 100 km for artificial satellites, is 
emphasized for future developments in the physics of the earth’s 
upper atmosphere. 
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THE POLAR UPPER ATMOSPHERE AND AURORA 


1. INTRODUCTION 

The upper atmosphere in the earth’s polar regions is the locale of many 
interesting phenomena and unsolved problems arising from the following ter- 
restrial facts: 

1 . The earth possesses a magnetic -dipole field whose axis is very close 
to its spin axis.* 

2. The earth’s spin axis is not normal to the plane of its orbit around the 
sun, but is inclined to it at 66.5 degrees. 

Because of fact 1, solar corpuscular radiations** have more effect in the 
polar regions than in the lower latitudes, causing the polar absorption (PC A) 
events, auroras, and related phenomena in the polar upper atmosphere. 

Because of fact 2, the polar regions have prolonged sunlit and dark periods, 
resulting in large seasonal variations. 

It should be noted, however, that the largest energy input from solar radia- 
tion (the main source of dynamical activities of the earth’s atmosphere) occurs 
in tropical regions where the diurnal variation dominates. In this paper, the 
following subjects are discussed; first, the so-called winter anomalies in the 
ionosphere and mesosphere; second, characteristics of the aurora with its re- 
lated phenomena and some aeronomical problems in the polar atmosphere; and 
finally, the latest developments in the studies of noctilucent clouds. 

The "polar region" discussed in this paper is not necessarily confined to 
the arctic and antarctic zones; in some cases it includes geographic latitudes 
below 60 degrees. Geomagnetic poles (especially the south pole) are displaced 


*The dipole center is displaced approximately 340 km from the earth’s center toward the direction 
6.5° N, 161 °E and moving northwesterly at present. 

**These should be called “solar plasma," except the high energetic solar protons which are called 
"solar cosmic-rays.” The plasma holds neutrality as a whole by an equal amount of proton and 
electrons. (Solar plasma has another name, “solar-wind.”) 



by more than 10 degrees apart from the earth's axis (geographic poles).* This 
has great longitudinal effect on geomagnetical events. 


2. WINTER ANOMALIES 

If the solar ultraviolet-radiation is regarded as the sole source of ioniza- 
tion and heating of the upper atmosphere, there are several phenomena which 
cannot be explained by this concept. These are called "anomalies." (Ratcliffe, 
1960). 

2.1 F2 Winter Anomaly 

The maximum radio frequency reflected by the earth's ionosphere is the 
critical frequency of the F2-layer,f 0 F2. This corresponds to the plasma reso- 
nance frequency at the F2-peak which is given by (Appendix A) 

f ^ 9 • 10 3 VlT, (2.1) 

m m ' 


where N m is the number of electrons per cc at the F2-peak. 

If one ignores the dynamical motion of air in the F-region, the daytime N m 
can be obtained by the equilibrium between the production rate due to air photo- 
ionization and the loss rate due to collisions (attachment and recombination). It 
can be easily shown that (Appendix B) 

_L j. 

N m “(cosX) 2 , or f 0 F2 co (cos X) 2 - / 2 . 2 ) 


where X is the solar zenith angle. 


*ln general; two types of poles can be defined: 

(1) Dip-pole (magnetic pole), the point on the earth’s surface where the magnetic field is vertical; 
(In 1965, their locations are 75.5°N, 101. 0°W and 65.5°S, 140. 3°E, respectively.) 

(2) Geomagnetic pole, which corresponds to the intersection of the centered dipole axis with the 
earth's surface (i.e. 75.5°N, 101. 0°W for north-pole in 1965). For global-scale problems, the 
geomagnetic pole is more important than the magnetic pole (the latter is greatly influenced by 
local magnetism in the earth’s crust). 
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This indicates that the maximum values of f 0 F2 (orN m ) and the amplitude 
of its diurnal variation should be larger in summer than in winter. Observed 
results, however, disagree with this conclusion (Figures 1 and 2). 

Note that the maximum values of f 0 F2 and its diurnal variation occur in the 
auroral zone (around 70°N or S) and decrease again in the polar cap region (Fig- 
ure 1) and both maxima (noon value and amplitude) increase with sunspot activity 
(Figure 2). These anomalies show further complications, such as universal- 
time (UT) dependence and hemispheric asymmetry (higher electron density in 
the southern hemisphere). 

These anomalies may be due to dynamical effects such as the following (for 
which there are several arguments): 

1. Trans-hemispheric meridional circulation (F. S. Johnson, 1960). 

2. Horizontal wind shear and geomagnetic control (Hill, 1960). 

3. Ionospheric current (Sato and Rourke, 1964). 

4. Ionization by solar-charged particles. 

Effect 1 brings electrons in the summer hemisphere (photo-electrons) by 
direct flow to the winter hemisphere, where the contraction of large-scale sub- 
sidence further enhances the relative electron density. This explains the higher 
F-region electron density in the polar winter, but not its large amplitude of 
variation. 

Effect 2 can be seen from Figure 3. This effect is combined also with Ef- 
fect 3 (Duncan, 1962). Apart from dip-pole, the magnetic lines of force are 
slightly oblique. If two winds flow in opposite directions, electrons accumulate 
between the two layers of counter-wind (wind shear). This is because electron 
motion is restricted along the magnetic line, while dominant neutral particles 
move with horizontal wind. Since winds and currents depend on the local time 
and season, the above explanations have been accepted so far but have yet to be 
proved by direct observations of wind or current. Note in Figure 3 that the posi- 
tion of the f 0 F2 peak should fit in these two layers and that, if wind directions 
are reversed, the peak between those layers is reduced rather than enhanced. 
According to the top-side sounder (Alouette 1, 2, Explorer 22) in the F-region, the 
winter anomaly is limited to the region below the F2-peak. Above the peak, the 
electron density is nearly two orders of magnitude higher in the summer polar 
ionosphere. This can be explained by higher temperature, and correspondingly 
larger scale heights, above the F2 peak at the sun-lit (summer) pole. 
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2.2 Winter Anomaly in the D-Region 

In the lower ionosphere, electrons are not only forced to oscillate by the 
radio-wave electric field, but also collide with neutral air molecules in the dense 
atmosphere, dissipating radio wave energies into the thermal energies of upper- 
air molecules. In other words, ionization in the lower ionosphere, below 90 km, 
causes radio-wave absorption. 

Absorptions in the D-region are associated with several other geophysical 
phenomena, which are discussed in the later sections, and have been observed by 
many different techniques such as (1) continuous measurements of cosmic noise, 
(2) forward scattering (oblique long-distance communications) , and (3) observa- 
tions by "ionosonde" — the ground-based vertical sounding of the ionosphere, 
giving f in , the ratio of multiple-echo intensities, etc. 

Based on ionosonde data for 1942-1949 at Slough, England, Appleton and 
Piggott (1954) reported the following findings: 

1. Maximum absorption occurs in winter. 

2. The amplitude of the diurnal variation is larger in winter than in summer. 

3. Maximum absorptions and amplitudes of diurnal variations both increase 
with increasing sunspot number. 

These can be seen from Figure 4 which plots the average monthly noon values of 
absorption against inverse Chapman function of solar zenith angle X. (The in- 
verse Chapman function is well approximated by cos X except when X exceeds 
80 degrees.) 

Two lines in Figure 4 correspond to the two values of the solar activities, 
i.e., R = 200 and 100, respectively, where R is the average annual sunspot num- 
ber. Figure 5 also shows these anomalies; monthly median values of absorption 
(in decibel) for conditions of constant solar zenith angle cos X = 0.25 (i.e. X = 75° 
approx.) at Slough are shown for the period 1942-1949.* Note that the constant 
solar zenith angle, cos X = 0.25, corresponds to the noon value in winter (at 
Slough) but to early morning and late evening in summer. 

The winter anomaly absorption was first noticed by Appleton and Piggott 
(1954). Since then, many statistical studies of this phenomenon have been made 


*The values of absorption in decibels are reduced to an equivalent frequency of 1 Mc/s (MHz) by 
an empirical formula. 
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by different workers (Dieminger, 1952; Beynon and Davies, 1955; Thomas, 1902), 
and the following additional features have been revealed: 

1. Days of high absorption in winter occur in groups, which show no cor- 
relation with geomagnetic activity. 

2. Occurrence of high absorption indicates a maximum near 55 degrees 
geographic latitude (over Europe). 

3. The winter anomaly is observed over an area of at least 10 6 km 2 
simultaneously. 

There are essentially two different interpretations for this anomaly: (1) the 
precipitation of charged particles; and (2) meteorological effects. The correct 
interpretation will be found in the future by direct observations (rocket- 
soundings and satellites). Meanwhile, the outlines of these two ideas seem worth 
discussing. 


2.2.1 Anomaly attributable to meteorological effects 

Association of the absorption enhancement in the D-region with increasing 
temperature in the stratosphere has been reported by many workers (Bossolasco 
and Elena, 1963, 1966; Gregory, 1961, 1965, 1966; Shapley and Beynon, 1965). As 
an example, Figure 6 shows the non-deviative absorption index A [db(MHz), Ap- 
pendix B] and temperature of the 10-mb level around Central Europe for two 
winter periods: 1958-1959 and 1960-1961 (Bossalasco and Elena 1963). 

It should be noticed that thermal expansion of the stratosphere does not 
compress its upper layer but rather reduces upper layer densities including 
ions in the D-region, causing outflows from the region where the upward ex- 
pansion has started. Moreover, it has been known that the warmings and the 
absorption enhancements are not always simultaneous; often there are 2 or 3 
days' lag either way. Until we know the cause of stratospheric warming, we 
cannot distinguish cause and effect (Kenecht, 1965; Belrose, 1963-1965). 

As shown by Nicolet (1965), the warming of the D-region can produce the 
enhancement of nitric oxide (NO) density. Because of its relatively low ioniza- 
tion potential, NO is regarded as a main ion constituent in the D-region.* 


*The threshold ionization potentials of N 2 , 0 2 and NO are 15.6, 12.2 and 9.2 ev, respectively; 
corresponding threshold wavelengths for photoionization are 820, 1026, and 1345 A, respectively. 
The only species that can be ionized by the solar Lyman alpha (Lya, 1215-7 A, 10.2 ev) is there- 
fore NO (Nicolet and Aikin, 1960). Recent rocket measurements by mass spectrometer indicate, 
however, that below 90-km level N0 + is a minor ion and major ions (over 70 percent) consist of 
large-mass species such as 37 + , 45 + up to 55 + , etc. (Narcisi and Bailey, 1965) 
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From the direct rocket observation that indicates the increase of nitric oxy- 
gen concentration below the 100- to 75- km region (Barth, 1964), Sechrist (1967) 
recently developed a theory to explain the D-region winter anomaly by the well- 
known warm winter-time mesosphere and by occasional temperature inversion 
in the winter mesosphere (Jones, et al, 1959). 


2.2.2 Effects of particle precipitation 

The first paper on the D-region winter anomaly presented by Appleton and 
Piggott (1954) is based on the comparisons of attenuation of a radio wave at con- 
stant solar altitude at different seasons (see Figure 5). As pointed out already, 
therefore, the attenuations are compared at different local times — noon in winter 
but near sunset or sunrise in summer. As shown in Figure 7, the precipitation 
of electrons is maximum around noon-time (L = 4 corresponds to geomagnetic 
latitude 60 degrees).* (The scattered points in Figure 7 must be examined with 
an eye on the "median" or trend line.) 

From these findings, Maehlum (1967) recently concluded that the total rate 
of ionization in the D-region, q, is the sum of irregular ionization caused by 
charged particles and regular ionization, q R . The latter is due to the solar 
ultraviolet-radiation and depends on the zenith angle of the sun, X, and the 
former depends on the local time, t, i.e., 

q = q R (X' h ) + q r (th) 
where h indicates the altitude. 


*A dipole line of force is described by the simple equation R = L cos^A, where R is the radial 
distance from the center of the dipole, and A is the geomagnetic latitude. The maximum value 
of R; i.e., radial distance in the equatorial plane, specifies the line of force and is called 
“L-value.” Three adiabatic invariants, constant along a given value of L, designate the motion 
of a trapped particle in the pure dipole field. The magnetic field in the presence of a magnetic 
dipole is 

B= — (1 + 3 sin 2 #, 

R 3 

where M is the dipole magnetic moment. The L-B coordinates are equivalent to the R-A coordi- 
nates. Because of the direct connection of L-B coordinates with the adiabatic invariants of par- 
ticle motions, L-B coordinates are the more commonly used to describe the data of trapped par- 
ticles in the earth’s magnetic field, even though the field is not that of a pure dipole (Mcllwatn, 
1961). 
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In summer q R predominates giving approximated q « cos X . In winter q r 
(which has strong local-time dependence and day-to-day variation) predominates. 

Appleton and Piggott (1954) concluded that strong absorption phenomena 
could not have been attributed to particle precipitation because they have no 
correlation with high geomagnetic activity, K p . However, according to Maehlum 
(1967), the large positive correlation between the flux of particle precipitation 
and K p does not necessarily conflict with a poor correlation between high absorp- 
tion and particle precipitation: to increase the absorption only needs a flux of 
high-energy electrons or hard-spectrum electrons. Since most high-flux pre- 
cipitations consist of low-energy electrons, the precipitation of high-energy 
electrons is independent of K p . An unsolved problem is, therefore, the relative 
importance of q R and q r . 


2.3 The Mesospheric Temperature in Winter 

In winter, that is, the polar night, the stratosphere is very cold; but, the 
high- latitude mesosphere (above 70-km altitude) is warmer than in summer. 

From studies of anomalous sound propagations, meteor trails, and especially 
the wind system above 60 km, this warm polar-night mesosphere was known 
before IGY investigations began. 

From intensive investigations of general atmospheric circulation, Kellogg 
and Schilling (1951) suggested that the adiabatic heating of polar -night meso- 
sphere is due to large-scale subsidence in the polar region. This warm meso- 
sphere in the polar night was confirmed by the direct rocket sounding during IGY 
(Stroud, Nordberg, Bandeen, Bartman, and Titus, 1959). In 1961, Kellogg intro- 
duced a new heating process, called "chemical heating" by the recombination of 
dissociated oxygen molecules. This drastically reduced the required velocity of 
subsidence which had been more than 1 km per day to compensate for cooling 
from infrared emission by carbon dioxide and water vapor. As pointed out in 
the previous section, this warm polar mesosphere in the winter hemisphere 
supports Barth's interpretation of the high mesospheric radio-wave absorption 
because of the high concentration of nitric oxide in warm air. 

Finally, the radiative heat source in the mesosphere has been shown to be 
the summer polar mesosphere; the heat source in the troposphere, is always in 
the equatorial region and the corresponding heat sinks are near the polar tropo- 
pause and over the winter pole in the mesosphere (Murgatroyd and Singleton, 
1961). It is, therefore, dynamical processes, i.e. , convection from the summer 
mesosphere to the subsiding polar winter mesosphere, that causes the anomalous 
warm mesosphere at the winter pole and the cold mesopause at the summer pole. 
This dynamical process can be regarded as the thermally driven annual oscillation 
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of the earth's atmosphere which has been also discussed theoretically by 
Sawada and Matsushima (1964). According to recent rocket soundings, the 
vertical temperature distribution in the warm polar mesosphere indicates wide- 
range daily fluctuation, while those of the cold summer mesopause is steady for 
many days (Theon, Nordberg and Smith, 1967). The cold summer mesopause is 
one of the important factors in the occurrence of noctilucent clouds (which are 
discussed in the last chapter of this paper). 


3. AURORAE AND RELATED PHENOMENA 

The aurora borealis, one of the most beautiful and fascinating phenomena in 
the polar upper atmosphere, has been observed and described by many scholars 
in northern Europe since the early 17th century. The complexity of this phenom- 
enon in connection with solar activities and the earth's magnetic variations has 
been discussed in various books such as those by Birkeland (1896, 1908, 1913), 
Alfven (1950), Harang (1951) and Stormer (1955). The discussions are based on 
visual observations from the ground, laboratory experiments, and on mathemati- 
cal calculations. Since World War n, research has been greatly advanced by 
direct rocket soundings and satellite observations. A detailed summary of pre- 
IGY auroral research is contained in a book by Chamberlain (1961). Further 
developments in auroral physics made by the use of world-wide records of all- 
sky cameras, magnetograms, and rocket and satellite observations are described 
in an extensive article in the Encyclopedia of Physics (Akasofu, Chapman, and 
Meinel, 1966). The following sections of this paper will discuss some important 
features of this very complicated phenomenon. For further details refer to the 
above sources. 


3.1 Temporal and Spatial Variations 

At about 60 degrees geomagnetic latitude* (just outside the so-called auroral 
zone**) in the evening of a geomagnetically disturbed day, a faint white greenish 
glow is visible over the northern horizon around 9 or 10 p.m. (local time). Grad- 
ually this glow spreads toward the south, covering the whole sky, which is now 


*To avoid redundant arguments, the description in the following sections is limited only for the 
northern hemisphere unless stated specifically. From symmetry, similar arguments should be 
applied to the southern hemisphere. 

**There is a maximum number of aurorae along 67 degrees of geomagnetic latitude. The number 
of aurora appearances decreases rapidly above 70 and gradually below 65 degrees of latitude. 
These annular zones (north and south) around geomagnetic poles, in which auroral activity is 
most frequent and strong, are called “auroral zones." As will be discussed later, “auroral oval" 
is more effective in describing individual (temporal and spatial) auroral activity. “Auroral zone 
corresponds to the trace of the midnight portion of auroral oval and can be defined statistically. 
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seen to consist of arcs stretching east and west. The number of arcs increases, 
drifting towards the west. If the geomagnetic disturbance is not great, the auroral 
display fades out into the western sky. On the other hand, if a geomagnetic storm 
has started (which can be easily detected by checking a rapid-run magnetogram), 
one of the arcs will brighten suddenly around the midnight sector in the auroral 
zone. Then the fantastic auroral displays begin, spreading over the whole sky. 

This most exciting phase of auroral activity is called "auroral break-up." Sup- 
pose that we were in the western side of this break-up, i.e., in the evening portion 
of the auroral zone; then we would see a sudden brightening in the eastern sky, 
and the arcs swiftly approaching. This fast movement of the arcs toward the west 
is called "westward surge." On the other hand, in the eastern side of the auroral 
break-up (or in the morning portion of the auroral zone), all arcs are spread in 
patchy forms, pulsating and drifting toward the east. This is called "eastward 
surge." The eastward drift is faster than the westward drift i.e. that of eastward 
surge is of the order of 10 km per sec; while the westward of the order of 1 km/ 
sec. The westward surge moves along pre-existing arcs; in the eastward surge, 
the auroral system disintegrates into irregular folds and patches. The brightness 
of these patches fluctuates rapidly and eventually fades intothe early morning twi- 
light. These variations differ case by case; sometimes cycles repeat twice or 
more in a day; i.e., an explosive expanding phase in which auroral activity spreads 
simultaneously toward lower and higher latitudes , and a recovery phase in which 
the northern extended edge retreats slowly and only quiet arcs remain. Figure 8 
shows the temporal variation of auroral activity observed from a single station as 
a typical example (Heppner, 1958). Note that the same displays are seen differ- 
ently from different locations. For example, arcs far from the observing station 
are seen as a glow because their lower borders are below the horizon. If the 
arcs are overhead, they appear fan-shaped, radiating from a point that corre- 
sponds to the geomagnetic zenith, and form a "corona." The morphologies indi- 
cated in Figure 8 should be read, bearing in mind the effect of location on the 
aspect of an aurora. By combining several pictures taken by all-sky cameras 
simultaneously along the Arctic circle (Canada, Alaska, and Siberia stations), it 
is found that the auroral arc forms a very large loop along a line called "auroral 
oval." Inside the oval, there are the so-called "polar cap aurorae," in which arcs 
and bands form fairly straight lines from noon to midnight meridian across the 
geomagnetic pole. These features can be seen from Figure 9, supplied by 
Akasofu (1964), which indicates the spatial distributions of auroral arcs and 
their temporal variations and directions of motion. Note that the widths of these 
arcs are only of the order of 200m,* and that parts of the arcs contract into small 
circles that are especially bright; these are called "rays." 


*The gyroradius of 20-kev electrons at these altitudes is of the order of 20m. Because of colli- 
sions with neutral molecules, however, horizontal spread should be larger than this in the earth's 
magnetic field. 
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In Figure 9, concentric circles indicate geomagnetic latitudes. Parts (a) 
and (b) correspond to a quiet state; part (c) shows the onset of auroral break-up; 
part (d) shows the peak of auroral activity, parts (e) and (f) indicate decaying 
phases. Directions of auroral movement are shown by arrows and the position 
of the sun (noon) is the top of each figure. 


3.2 Geomagnetic Storms and Aurorae 

The relation between auroral activities and geomagnetic variations has been 
investigated by many workers (Heppner, 1958, for example). Relations between 
the so-called ’’geomagnetic bays"* have been examined extensively by Akasofu 
(1965), using world-wide data from a magnetogram and an all-sky camera. Con- 
ventionally, the geomagnetic bays are explained by the electric -current system 
and its time variations, which are assumed to be flowing in a maximum conduc- 
tive layer in the lower ionosphere around 100 km (i.e., lower E-region) level. 
Except for the daytime equatorial ionosphere, these electric currents are not 
directly observed, and should be regarded as an equivalent current system, be- 
cause these currents might not necessarily exist solely in a thin layer as is con- 
ventionally assumed. The ionosphere current produces a positive bay (horizontal 
magnetic field increase) if the current flows eastward, and a negative bay if the 
current flows westward. Akasofu (1965) used these observational findings to 
propose a new scheme of the equivalent current system in the polar (or "auroral") 
substorm (also called "auroral substorm"). This is shown in the Figure 10(b) in 
comparison with the previous model (a). In part (b) the westward current (auroral 
electrojet) flows along the auroral oval. 

These two models of the equivalent current system stand for a maximum 
phase of polar substorm. The cycle of auroral substorm can be repeated several 
times in a day during a period of high geomagnetic activity, which consists of a 
sudden brightening of arcs and explosive expansion of the active area, within 10 
minutes or so, and the slow decaying phase that lasts for one or two hours. 
Breakup (in which all arcs or bands disintegrate into very bright patches) occurs 
only around midnight. In the previous model (Chapman and Bartels, 1940; 
Fukushima, 1953; Vestine, 1940 and Nagata and Fukushima, 1952), shown in Fig- 
ure (10a), it is assumed that the equivalent current (in either direction) flows 
along a geomagnetic latitude (in the auroral zone) and the equivalent west- 
ward current flowing in the nightside is more intense than the eastward cur- 
rent in the late morning (daytime) side. On the other hand, in Akasofu's new 


*ln a record of magnetometer (horizontal component), a bay-type variation appears for a short period 
of the order of hours. This is called the “positive bay" if variation shown a maximum and the 
“negative boy" if variation shows a minimum. 
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system the equivalent current flows westward along the auroral oval whose 
night side position is close to the auroral zone, but the daytime position is 
shifted into the polar cap. At point A, both current systems produce a strong 
positive bay, and at B a strong negative in both systems. The new system, 
however, produces a weak positive bay at C, which is in the evening sector out- 
side of the auroral zone. The previous system produces a weak negative bay at 
point C. 

The above-mentioned current system analysis can be summarized as fol- 
lows: (Chapman and Bartels, 1940, p. 301). 

The storm -time geomagnetic variation, D, can be expressed by 

D = Dst + DS + Di , (3.1) 


where Dst is the universal storm-time variation attributed to the (westward) 
ring current in the earth’s magnetosphere, 

DS is the "disturbance longitudinal inequality," that is, the function of 
geomagnetic longitude, 4> , and latitude, k , consisting of diurnal and 
semidiurnal components, and 

Di is the irregular part of D, dominant in the auroral zone, and is called 
"polar magnetic substorm" (Akasofu, 1966), 

Akasofu's new current system corresponds, therefore, to 


D - Di , 


(3.2) 


whereas the previous conventional current system corresponds to 


and 


D - Ds t + DS 

(3.3) 

Dst = C 0 (k) 

(3.4) 

(k) sin (n<p + 2n) n > 1 


(k) sin (ncfi + 2j) 

(3.5) 


Note that the pattern corresponding to DS (Equation 3.5) is fixed with respect 
to the sun and shifts westward with the rotation of the earth. As pointed out by 
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Akasofu (1966), however, in the polar region a current system resembling SD 
appears intermittent with a lifetime of the order of 1 to 3 hours. 


3.3 Spectral and Aeronomical Aspects 

The aurora brightness described in the previous sections (shown in pictures 
of arcs taken by an all-sky camera) is mostly due to emissions of the forbidden 
line (5577 A), a transition from metastable oxygen atoms, ( X S 0 - *D 2 ) which ap- 
pears as a greenish white glow. The vertical distribution of this brightness is 
approximately a Chapman function with a peak around the 100-km level; i.e., its 
intensity decreases nearly exponentially with altitude above the maximum and 
decreases very rapidly below it (the layer of maximum brightness). There are 
red aurorae that have been classified into two types: one, "Type B," appears at 
a very low altitude mostly at the bottom of bright arcs, as a purplish red border. 
The emission of type A aurorae consists of two lines from a neutral oxygen atom 
01 (6300 A, 1 D 2 - 3 P 2 and 6320 A, 1 D 2 - 3 P, ) (see Figure 11) and of the first 
negative band* from an ionized nitrogen molecule, N* (3914 A, B 2 2 + - X 2 2 + ). 

Since these red lines from the neutral oxygen atom are due to the transition 
from a very long-life metastable (110-second) state, 1 D 2 , to the ground states 
3 P 2 and 3 Pi , heights of emission are limited above 200 km where deactivation** 
by collision is negligible. On the other hand, Type B emission is mostly due to the 
first positive groups of neutral nitrogen molecules (5854 A ~ 10440 A, B 3 n g _ 

A 3 2 u ) (Figure 12). The units of brightness, energetic relation, and energy dissi- 
pations are discussed in Appendix C. Figure 11 indicates main transitions in 
oxygen and nitrogen atoms for neutral, [OI], INI] and singly ionized ions [Oil], 

(Nil] . Figure 12 indicates main transitions for molecular oxygen (Qj , 0 2 ) and 
nitrogen (N 2 , N 2 ), respectively (Davis, 1965). 

Since the height of maximum intensity of Type B aurorae is around 80 to 
90 km, energies of the incident primary electrons that are known as the source 
of this type of aurora should be higher than 200 kev. Although the brightest 


*The names of positive and negative bands (or groups) originate in cathode-ray or glow-tube 
experiments. Spectra of the glowing columns near the anode and cathode are called positive 
and negative groups, respectively. These correspond to the emissions from neutral and singly 
ionized ions, respectively. 

**Transition to other states without photo-emission is called “deactivation,* where the energies 
are balanced by colliding particles (molecules, atoms, or electrons) instead of photons. 
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aurora is due mainly to the forbidden 01 (5577 A, line, there are many permitted 
lines from neutral oxygen, neutral nitrogen, and ionized nitrogen atoms. None of 
these lines is however, stronger than the forbidden line. 

Note that the above-mentioned visible emissions are all due to auroral elec- 
trons impinging from outer space on the earth's polar upper atmosphere. There 
are also visible light emissions due to auroral protons; these emissions can be 
classified in two categories: (1) hydrogen emission and (2) PCA or polar-glow 
aurorae. The hydrogen emissions observable from the ground are essentially 
Balmer lines Ha (6562.8 A), H (4861.3 A) and Hy (4340.5 A ; , which appear in 
the quiet homogeneous arcs during the early stages of auroral activity at night. 
Although the total intensity of hydrogen emission is far weaker than that of the 
electron-induced aurorae, the discovery of Doppler shift in the H a profile by 
Meinel (1951) brought great progress in auroral physics. Doppler shifts of 
Balmer lines toward the geomagnetic zenith direction are toward the violet 
side and asymmetric, whereas it is fairly symmetric in the geomagnetic hori- 
zontal direction. The asymmetric Doppler shift gives the velocity of impinging 
protons, which is of the order of 3000 to 4000 km per sec. This is consistent 
with those energies estimated by the range-energy relation of protons in air, i.e., 
around 300 kev to penetrate to the depth around 100 km. 

The following two reactions can be regarded as the sources of the excited 
hydrogen atom, H*, which emits Balmer lines: 


+ H - H* + H + , 

(3.1) 

+ 0 - H* + 0 + , 

(3.2) 


where H + is the primary auroral proton entering the atmosphere. 

The neutral atomic oxygen, O, is more than a million times abundant than 
the neutral hydrogen atom, H, around the altitude 100 km, and since the ioniza- 
tion potential of oxygen atoms (13.62 ev) is very close to that of hydrogen (13.60 
ev), the charge exchange (resonance-type) of auroral protons with atmospheric 
oxygen atoms (Equation 3.2) is a main source of excited-state hydrogen, H*, the 
emitter of Balmer lines in aurorae. Above 500-km altitude, another charge- 
exchange reaction (Equation 3.1) becomes important, which has been known as a 
possible source of the ultraviolet Lyman series. The cross sections for reactions 
3.1 and 3.2 are very large below 10 kev (both of the order of 10~ 15 cm 2 ), but de- 
crease rapidly with energy above 10 kev (less than 10 -18 cm 2 at 200 kev). Emis- 
sion of hydrogen lines occurs, therefore, near the end of the range of incoming 
protons; i.e., the emission dominates below a 100-km level (Eather, 1967). 
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Note that these temporary neutralized auroral protons spread horizontally 
more than several hundred kilometers across the earth's magnetic field before 
being ionized again or stopped by colliding air molecules. Davidson (1965) has 
shown by means of Monte Carlo type calculations that these spatial spreads of 
auroral protons occur in the polar upper atmosphere. The results explain why 
auroral hydrogen emission does not have a fine structure such as arcs or bands 
produced by auroral electrons. Note that the hydrogen emission lines are not 
only in the visible region, but also in the ultraviolet (Lyman) and in the infrared 
(Paschen, Brackett, etc.). Because of the atmospheric ozone, the Lyman lines 
are not observable from the ground; they have been measured by rocket sounding. 

Another proton-produced aurora - the "polar glow aurora" — is a vast air 
glow occurring simultaneously with so-called "polar-cap absorption" (PCA).* 
According to Sanford (1962) , the maximum brightness is comparable to that of 
the Milky Way, although strong emissions are essentially the same as the type B 
(purplish red) aurora.** This similarity to Type B aurora is predictable, be- 
cause the polar-cap glow is due to the direct excitation and ionization of molecular 
nitrogen and atomic oxygen by the solar protons which penetrate the atmosphere 
below the 100-km level, (proton energy is 5 to 100 Mev). 

Figure 13 shows the example of the polar glow aurora reported by Sanford 
(1962), indicating the time variations of auroral glow, with cosmic-noise absorp- 
tion (riometer data at Thule) , K-index at Scott station and Kp-index. 


3.4 Auroral Infrasonics 

The large-scale auroral displays produce sound waves by disturbing the 
earth's tenuous upper air.*** These waves have frequencies below audible range 


*After the eruption of flares on the sun’s surface, mostly within a few hours, signals of the 
high-latitude-path radio wave disappear. This is called “polar cap blackout" and has been 
explained by D-region ionization enhanced by high-energy solar protons (5 to 100 Mev). By 
means of a “riometer* (a relative ionosphere opacity meter that measures relative-intensity 
variations of cosmic-noise around 30 MHz), details of D-region ion density variation have been 
studied and named “Polar Cap Absorption." 

**During intense auroral activity, auroral bands appear across the polar cap, as mentioned before. 
These bands are produced by electrons, whereas the polar cap glow (a weak glow covering whole 
polar cap area almostuniformly) is produced by protons. Thus, both are clearly distinguishable. 

***Audible sounds during auroral activity have been reported by many observers. These seem rather 
local phenomena, if it is real, probably due to the electro-static discharge near the ground through 
the dry, cold air. (Chapman and Bartels, 1940, p. 466) 
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and are called "infrasonic waves" or "infrasonics." Since there are many kinds 
of pressure waves and disturbances in the atmosphere, especially near the ground, 
these subaudible sounds with intensity of the order of several dyne/cm 2 * should 
be listened to with the aid of the following device: 

In order to eliminate non-periodic local disturbances such as the pressure 
fluctuation due to wind, a long pipe with many holes spaced at random intervals 
is attached to each detector's input (microphone). Local pressure fluctuations 
whose scale is less than the length of the pipe can be cancelled out; only large- 
scale waves can be detected. In order to find the direction of propagation of the 
traveling pressure wave, more than three of these detectors are installed approx- 
imately 10 kilometers apart. The differences between the times when the wave 
reaches the detectors indicate the normal direction of the traveling wave front. 

Figure 14(a) indicates the diurnal variations of arrival direction of infra- 
sonic waves observed at the National Bureau of Standards (NBS) stations** in 
Washington during magnetic storms, as measured by the above procedures. As 
described in section 3.1, auroral activity shifts along the auroral oval within a 
day, and the most active part, (i.e., the auroral breakup occurring around the 
midnight sector) moves essentially along the auroral zone clockwise around the 
geomagnetic pole. Three cases are shown in Figure 14(b), corresponding to 
three different local times in Washington*** (21, 0, and 3 o'clock, respectively); 
the hatched area in each case indicates the position of the most active part of the 
aurora. The location of Washington is also shown by an asterisk in front of an 
arrow. The time-variation of the arrival direction of the traveling pressure 
wave (infrasonics) during a magnetic storm observed in Washington agrees with 
the assumption that these waves originate in the active part of the aurora. Since 
the pulsating aurora appears just after the auroral breakup (Figure 8) it was 
thought that the source of the waves might be the pulsating aurorae caused by 
fluctuating ionospheric heating from the precipitations of auroral electrons 
occurring at the most active phase of aurorae (Maeda and Watanabe, 1964). 

According to the recent observations performed inside the auroral zone, 
the sources of these infrasonics are not pulsating aurorae but the supersonic 
horizontal movements of auroral arcs (Wilson and Nichparenko, 1967). As de- 
scribed in the section 3.1, auroral arcs move rapidly westward in the evening 
sector (i.e., the west side of the breakup) and shift eastward in the morning 
sector (i.e., the east side of the breakup). This is also shown in Figure 9, where 


*0f the order of one millionth of atmospheric pressure on the ground. 

**The details of equipments and locations of stations are reported by Chrzanowski et al. (1961) 

***75° WMT. 
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the arcs shaped like a capital omega spread both westward and eastward around 
the midnight section, (i.e. , the breakup portion), constituting the westward and 
eastward surge. As mentioned in section 3.1, the eastward surge is generally 
faster than the westward; both are generally supersonic. Thus, shock waves can 
be formed and propagated in the normal direction of these arc-fronts. Figures 
15 and 16 show an example of the detailed correspondence between the appearance 
of auroral infrasonic waves and the supersonic auroral surge. Figure 15 shows 
the temporal change of the position of the auroral arc, projected on the map at 
1-minute intervals (the time indicated at each arc-position is Alaskan Standard 
Time (AST) — the local time at 150°W). Figure 16 shows infrasonic records 
corresponding to this movement of the auroral arc; this is the superposed record 
obtained at three stations near College, Alaska. The time marks in the horizon- 
tal scale are 1-minute intervals; numbers are in UT (UT = AST + 10). A vertical 
bar at the right indicates the scale of 5 dyne per cm 2 ; 0 and V are the azimuth 
of wave normal (measured from north clockwise in degrees) and the horizontal 
velocity (in m per sec), respectively. Note that 0 and V are determined by 
shifting the records with respect to each other until the clear wave fronts in 
each record fit as shown in Figure 16, where time differences are omitted be- 
cause they are of the order of several seconds. Since V is the horizontal pro- 
jection of wavefront velocity (i.e. V~ c sec 9 , where Q is the angle between 
wavefront and the horizontal and c is the velocity of sound) , the altitude of the 
source can be deduced from these observations. The results agree with the 
assumption that the altitude of the source is around 100 km (the altitude of the 
bottom auroral arcs). Note also that the horizontal motion of auroral arcs and 
bands does not mean the motion of a real mass of glowing air, but a visual effect 
corresponding to the horizontal shift of frontal area in the upper atmosphere 
when bombarded by the auroral electrons. In this respect, the essential mecha- 
nism for pressure-wave generation is auroral heating (Maeda and Watanabe, 

1964; Cole, 1965) rather than the ion- motions and their collision with central air 
caused by the electrodynamic drift of electric current in the lower ionosphere 
such as the "Pedersen drift" (Piddington, 1964). Since the origin of auroral 
electrons is not known, the mechanism of supersonic horizontal drift of auroral 
arcs is not well understood. Another unanswered question is causes of the 20- 
to 80-second periodicity in the pressure wave, if rapid horizontal movement of 
arc (i.e., westward and eastward surge) causes the auroral shockwave that 
decays into an infrasonic wave in the lower atmosphere. 


4. NOCTILUCENT CLOUDS 

Another interesting and beautiful phenomenon in the polar atmosphere is the 
occasional appearance of tenuous, silvery blue clouds in the dark rim of pro- 
longed twilight in the summer sky. Since the unusual height of these clouds was 
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disclosed by Jesse (1896) these mysterious clouds, now called "noctilucent 
clouds,"* have been studied by many geophysicists (Vestine, 1934; Stormer, 
1935; Paton, 1949, 1950, and 1954). Only recently, however, has the mystery 
of these clouds been unveiled by world-wide general observations and by direct 
rocket observations (Fogle, 1962, 1964, and 1965; Paton, 1964, 1965, Fogle, 
Chapman and Echols, 1965; Soberman, 1963; Witt et al., 1964). 


4.1 Characteristics 

4.1.1 Features indicating atmospheric motions 

One of the striking features of noctilucent clouds is their great altitude, 
which has been measured more than 500 times over the past 80 years at vari- 
ous places in the northern high latitudes, and is known to be nearly constant 
around 82 km. These clouds generally have interwoven billowy patterns, also 
wave motions, that have been used to determine the thickness of the clouds 
although they are very tenuous. According to Witt's measurements (1964) the 
clouds are about 1 km thick, or less; the vertical amplitude of wave motion is 
around 2 km and occasionally 5 km. Measurements of cloud drift velocity have 
been reported since 1885. The most reliable determinations of drift velocity 
and wave motion have been made by time-lapse photography and motion picture 
technique (Fogle, 1966). Fogle has summarized old measurement records and 
his own observations are as follows: 

1. Drift velocity is mostly of the order of 50 m per sec with occasional 
high winds up to 200 m per sec. 

2. Drift is generally southwest in the northern hemisphere; northwest in 
the southern hemisphere. 

3. Wavelength (horizontal) ranges from 5 to 100 km with a preponderance 
maximum around 10 km. 

4. Wave propagation velocity derived from measurements of wave-crest 
movements, i.e., band movements, is more than 200 m per sec relative 
to the drift of the cloud system. 


*According to Fogle (1966), the name “noctilucent cloud” derives from 0. Jesse, whose systematic 
observations of these clouds date from 1885; he named them “Leuchtenden Nachtwolken.” This 
should be distinguished from “nacreous clouds” or mother-of-pearl clouds which appear at an 
altitude of 15 to 35 km d uring the winter. 
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Note that the noctilucent clouds are observed mostly around midnight. The 
above-mentioned drift motion of clouds should therefore be regarded only as the 
wind at the polar summer mesopause around midnight. 

Haurwitz (1964) interprets short waves as interface waves and long waves 
as internal-gravity waves. Propagation of the internal gravity wave is not strictly 
horizontal and its oblique propagation tends toward vertical with decreasing ve- 
locity (both phase and group velocity in the horizontal direction approach zero) 
as the wavelength decreases (Hines, 1960; Eckard, 1960; Maeda and Young, 1966). 
On the other hand, for waves longer than 10 km to constitute the interface wave* 
requires a very large vertical wind shear (or large discontinuity in the vertical 
distribution of temperature). It seems most plausible, therefore, that short waves 
below 10 km appearing in the noctilucent cloud are the interface waves due to wind 
shear around mesopause, and the longer waves are the internal-gravity waves. 

To confirm these interpretations requires further direct measurements of wind 
and of temperature distribution around mesopause. 


4.1.2 Optical features 

Since the discovery of their great height, which is far above a balloon's ceil- 
ing altitudes, the formation of visible clouds in such a high, tenuous atmosphere 
has been a mystery; many meteorologists have wondered whether their nuclei are 
dust or ice. 

Spectroscopic observations have shown that the spectra of noctilucent clouds is 
continuous with no emission lines , with stronger intensity in the blue than in the red and 
have several Fraunhofer lines. Although these observations indicate that the size of 
cloud particles is comparable with the wavelength of visible light, more definite 
estimation of particle size has been made by polarization-measurements of scat- 
tered light from the cloud. In Figure 17 the percentage polarization, P, of blue 
light - 4900 A observed at several occasions from the ground is plotted against 
the scattering angle, ^ (the angle between the directions to the scattering point 
in the cloud and to the sun — Witt, 1960, a and b). Two kinds of theoretical curves 
are also plotted for comparison; one is based on Rayleigh-scattering which is 
valid when the size of scattering particle is smaller than the wavelength of scat- 
tered light (shown by full lines); the other is based on Mie's theory, which depends 
not only on particle size but also on particle refractive index (shown by dashed 
lines). With (1960 b) made a similar calculation, assuming that the particle is not 
coated by a simple dielectric substance . The results of the latter are quite different 


*This is also called the “Helmholtz wave," which can be seen in billowy clouds in the troposphere 
(mackerel sky). 
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from observed values (Figure 17), especially at scattering angles exceeding 70 
degrees. Comparison of the observed values with Mie's theoretical values there- 
fore leads to the conclusion that the particles of noctilucent clouds are dusts of 
sizes around 0.1m, coated with ice, making a total radius of 0.2 to 0.5m. Finally, 
the most decisive observation, which settled these arguments was made from 
direct samplings of cloud-particles obtained by firing rockets into the clouds 
(Soberman, 1963; Witt et al., 1964). Two complete samples were obtained; one 
from visible noctilucent clouds, and the other from the same altitude but in the 
absence of clouds. Electron microscopic analyses* of these samples indicate 
that nickel and iron particles of the order of 0.1 M diameter are main constituents 
of the cloud-nuclei; they are coated by ice (possibly hoar-frost). This is con- 
firmed by the calcium film indicator of moisture used with sample collectors on 
rockets. These results agree, therefore, with the estimations based on polariza- 
tion measurement and its theoretical analysis. It is also found that particle con- 
centration is at least a thousand times higher in the presence of clouds than in 
their absence (roughly 10 particles per cc as against less than 0.01 per cc). The 
size and substances of the cloud nuclei suggest that these are of extra-terrestrial 
origin.** 


4.2 Theory of Cloud Formation 

The noctilucent clouds are observable only in the polar summer twilight 
when the sun is below the horizon with the solar depression angle a (Figure 18) 
between 6 and 16 degrees*** (Fogle and Haurwitz, 1966; Paton, 1964; Fogle, 1966). 


*Further measurements of characteristic X-ray spectra are made, exciting the nuclei by an 
electron beam of 10 to 30 kev; also, neutron activation analyses are applied, confirming that 
cloud nuclei are of iron-nickel-meteorite origin. 

**Details of rocket sampling and various types of analysis of the collected sample are reported 
in the six successive papers in Tellus Vol. 28, 1964, from page 84 to 117. Further discussion 
based on the size spectrum of dust nuclei is presented by Webb (1965, 1966 a) who concludes 
the dust is of meteorite origin. This conclusion agrees with the latest analysis of cosmic dusts 
and zodiacal light (Singer, 1967). 

***As can be seen from Figure 18, if a<6°, the background sky is too bright for the cloud to be 
seen from the ground P; if a> 16°, the sunlight reflected by the cloud is too weak to be seen. 
The earth’s shadow which forms a dark rim of the cloud is not that of the solid earth but of 
approximately the surface of 30 km in the atmosphere, as shown by H sc (Fogle, 1966). 
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The shaded areas in Figure 19(a), (b), and (c) - corresponding to geographic 
latitude, A. , of 40°, 60°, and 70°, respectively — show that the observable period 
of noctilucent clouds given by the above-mentioned geometric conditions varies 
markedly with k and the season. These figures indicate that the noctilucent 
clouds could be observed over a whole year at latitudes below about 60 degrees. 
However, Fogle (1966) has shown from his own surveys of IGY-network reports 
that noctilucent clouds occur most frequently at about k = 60°; negligibly below 
40° or above 70°.* On the other hand, analyzing the available data of upper air 
soundings, Murgatroyd (1P57) found the lowest temperature (at 80 km in summer 
high latitudes) to be lower than the radiatively predictable temperature under 
weak but continuous insolation in the polar summer. As shown in Figure 20, the 
cold- in- summer and warm-in-winter mesosphere has also been observed by 
rocket-grenade soundings at Fort Churchill ( \= 59°) (Stroud et al., 1960). 
According to the latest observations, the summer polar mesosphere is as cold 
as 135°K (Nordberg et al., 1965). By calculating global atmospheric circulation 
(which is consistent with observed wind systems in the upper atmosphere) Murga- 
troyd and Singleton (1961) have shown that there is a continuous ascending flow in 
the summer polar mesosphere, at 50 km to 80 km altitude, of the order of 1 cm 
per sec in the regions of noctilucent clouds. 

From his analysis of the mesospheric data obtained by the Meteorological 
Rocket Network (MRN), Webb (1965 and 1966) has also shown the global atmos- 
pheric circulation, indicating the coldest summer polar mesopause and ascending 
flow in the regions of noctilucent clouds. Based on this analysis, he presents a 
morphology of noctilucent clouds. His main points are as follows: 

1. Size and density distribution of the noctilucent cloud nuclei (10" 6 - to 
10" s - cm radius and 0.1 to 10 per cm 3 , respectively) is consistent 
with meteoroid influx (Whipple, 1963). 

2. Enhancement of meteoroid summer influx due to the forward-oriented 
northern hemisphere in the earth's orbit is another support for the 
meteorite origin of noctilucent cloud nuclei, (Fogle and Haurwitz, 1966). 
Note that these micrometeorites fall very fast until they hit the earth's 
dense atmosphere at altitudes called "meteor heights" (70-100 km - 
Roper, 1966) where meteors from outer space evaporate or disintegrate 
into micrometeorites. The terminal fall velocities of these micro- 
meteorites below these levels decrease rapidly to the order of several 
10 cm per sec. 


*He also shows that the maximum frequency of noctilucent cloud occurrence is earlier at lower 
latitudes; i.e., it occurs in early July for k& 50°, late July at 60° and early August for 

A.> 60 ° ) . 
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3. Viewed on a global scale, most water vapor diffuses into the stratosphere 
at the equatorial tropopause.* This moist stratosphere air is transported 
toward the pole by the poleward meridional circulation that predominates 
in the summer hemisphere at around 50 km** altitude. 

4. Because of oblique insolation at the higher latitude, the height of the region 
of effective heat input increases with latitude. Correspondingly, the strato- 
pause and stratospheric circulation shift upward from the tropical atmos- 
phere to the summer pole atmosphere. 

5. Because of continuous heating such as solar ultraviolet absorption by 
atmospheric ozone, there is a stable anticyc Ionic circulation in the sum- 
mer polar stratosphere. 

By making a superposition of the atmospheric diurnal tide on this anticyclonic 
circulation, Webb (1966) has shown a strong easterly wind in the night side of the 
summer polar stratosphere, and called it "stratospheric tidal jet." This is shown 
in Figure 21. The flow pattern in this figure, however, has no theoretical founda- 
tion. It should be noted that according to the theory of the thermally driven diurnal 
tide of the earth's atmosphere, amplitudes of tidal motions decrease in the high 
latitude although they increase with altitude (Lindzen, 1967). As pointed out in 
section 2.3, the cold-in-summer mesopause as well as the warm-in-winter 
mesosphere can be regarded as the result of the thermally driven global scale 
annual oscillation of the earth's atmosphere (Sawada and Matsushima, 1964). The 
anticyclonic and upward flow in the midnight section of the summer polar meso- 
sphere agree very well with the general southwestward movement of noctilucent 
clouds and the suspension of cloud-nuclei. 

Chapman and Kendall (1966) develop a more rigorous mathematical theory 
of the noctilucent clouds, starting from the following assumptions: (1) the cloud 
particles consist of water vapor condensed on dust nuclei; (2) the dusts are of 
meteoric origin, descending slowly by diffusion under gravity until they hit the 
"turbopause,” (the upper boundary of the turbosphere where turbulent motion 
predominates), then dispersing quickly toward the lower atmosphere; (3) water 
vapor has been brought up by convection to the mesopause; above the mesopause 
water vapor content decreased rapidly with altitude due to the dissociation by 


Thus, trie effectiveness of the cold trap 


at the tropopause and consequent dry stratosphere are 
questioned. The tropopause cold acts only to regulate the diffusive flow and thus increase the 
time required to achieve equilibrium in the moist lower stratosphere (Webb, 1965). 


*This height and flow velocity increase at higher latitudes. The average northward velocity is 

3 m per sec at K - 20°, but exceeds 10 m per sec at K = 60°. 
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solar ultraviolet radiation; (4) the turbopause occasionally descends from its 
usual height of 100 km to the mesopause (around 80 km), where brought-up 
water vapor condenses on the nuclei, forming the cloud particles that are seen 
by reflected sunlight. 

A mathematical discussion is outlined below. 

Since the sizes of falling dust-nuclei are less than the mean free path of air 
molecules (which is of the order of 1 cm at 80 km) Stokes' formula cannot be 
used to calculate their velocities of fall. Diffusion of dust particles whose radii 
are of the order of 10 - 6 cm in the air is discussed by calculating the "vertical 
diffusing flux." The diffusing flux consists of three parts: (1) ordinary diffusion, 
which tends to reduce the concentration gradient and make the composition uni- 
form; (2) gravity diffusion, by which the heavier element descends and the lighter 
one ascends; and (3) thermal diffusion, by which particles in the warmer region 
diffuse into the colder region, producing the highest concentration in the coldest 
part of the system. 

Since temperature increases with height above the mesopause, thermal dif- 
fusion tends to cause descent of dust particles. It is shown, however, that particle 
descent is due predominantly to gravity.* It is interesting to note that the Stokes- 
Cunningham equation (which is the Stokes' equation with a correction factor at- 
tributed to "slip of gas at the surface of the dust particle,"**) gives essentially 
the same value of fall speed as the gravity diffusion equation. *** 

Since the molecules of water vapor are lighter than average air molecules 
and since the atmospheric temperature decreases with altitudes below the meso- 
pause, both gravity diffusion and thermal diffusion cause the water vapor to 
ascend to the mesopause. For the mesosphere with a lapse rate of -5°C per km, 
it is shown that the gravity-diffusion rate of water vapor is of the order of 


*For example, in the thermosphere with temperature gradient 5°C per km^at rnesopause temper- 
ature 150° K, gravity-diffusion speed for dust particles of radius 2.10 cm and density 3 
(specific weight) is 18 cm per sec downward whereas the thermal-diffusion speed is only 
3.3 x 10 cm per sec downward. 

* * T h i s is because the molecular mean free path is comparable with or greater than the size of 
the dust particle. 

***Writing the radius and density of the dust particle as r and p, respectively, fall speed at the 
noctilucent cloud level is 3 x 10® r p cm per sec by the gravity diffusion and 2 x 10® x p cm 
per sec by the Stokes-Cunningham formula, respectively (Fogle, 1966). 
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2 x 10” 2 cm per sec, (i.e. , 17.3 m per day) and the thermal-diffusion rate of the 
order of 10~ 3 cm per sec. Because of the change of thermal gradient above the 
mesopause, upward diffusion of water vapor decreases. Furthermore, photo- 
dissociation by solar ultraviolet radiation increases with altitude. Therefore, 
as shown by Bates and Nicolet (1950), the concentration of water vapor decreases 
rapidly at altitudes between 80 km and 95 km, giving the scale height of the order 
of 2 km. 

When the turbopause hits the mesopause, the dust particles are dispersed by 
eddy diffusion. The coefficients of eddy diffusion in these altitudes are still the 
subject of speculation. The assumption that D E = 10D below mesopause and zero 
above mesopause, (where D E is the coefficient of eddy diffusion and D is the 
coefficient of molecular diffusion) implies the formation of a sharp ledge of dust 
at the mesopause. The results of computations done by Chapman and Kendall 
(1965) can be therefore summarized as follows: (1) because of the increasing 
photodissociation by solar ultraviolet light, the thickness of the upward- diffusing 
water-vapor layer is about 1.5 km; (2) because of eddy diffusion, dusts and nocti- 
lucent cloud particles are immediately dispersed below the mesopause, forming 
a sharp visible ledge of the clouds; (3) the upward-convecting velocity of 18 cm 
per sec to maintain the dusts at the mesopause can be reduced by hoar frost or 
flake-type water coatings of 10 -6 -cm; and (4) noctilucent clouds can appear only 
when the turbopause (whose usual height is around 100 km) is lowered to the 
mesopause at about 80 km. The mechanism of lowering the turbopause is left 
open. As pointed out by Webb (1966), lowering of the mesospheric temperature 
to saturate water vapor might be another important factor controlling the appear- 
ance of noctilucent clouds, it still seems necessary to measure temperature dis- 
tribution, turbulence, and the concentration of cloud-particles. 


4.3 Noctilucent Clouds and Aurorae 

By time-lapse photography, Fogle (1966) has shown not only the motion of 
noctilucent clouds but also the occasional appearance of aurorae and their effect 
on noctilucent clouds. The general effects are: (1) a reduction of the brightness 
and the extent of noctilucent clouds; (2) transformation of a well-ordered struc- 
ture into swirls or formless veils; and (in some cases) (3) disappearance of the 
noctilucent clouds. These observations indicate that the precipitation of auroral 
particles causes turbulence in the layer of noctilucent clouds and, in an acute 
case, might cause even the heating of the cloud level. As discussed in section 
3.4, rapid variations of auroral electron precipitation causes dynamic disturbances 
in the mesosphere. If these disturbances exceed a certain level, it is highly prob- 
able that the cloud particles around mesopause will be rapidly dispersed by eddy 
diffusion as shown by the computations of Chapman andKendall (1965). Subsequently, 
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the luminosity of the cloud diminishes. For auroral electrons to cause direct 
heating of the cloud layer, their energy should exceed the order of 100 kev. The 
precipitation of these high-energy electrons is not common in the usual bright 
aurora. To estimate the heating of the cloud layer by precipitating auroral elec- 
trons, it seems desirable to make direct measurements of atmospheric tempera- 
ture around the noctilucent cloud level. If temperature variation is negligible 
during the diminishing phase of noctilucent clouds following auroral activity, the 
coefficient of eddy diffusion around the noctilucent cloud level can be derived 
from its time variations. 


5. Summary 

Due to prolonged dark-night, twilight and daytime, there are peculiar annual 
variations imposed on the diurnal change in the polar atmosphere. Since the 
geographic poles are close to the geomagnetic poles, the effects of charged par- 
ticles emitted from the sun are far stronger in the polar upper atmosphere than 
at other latitudes. Most of these charged particles follow the magnetic lines of 
force around the earth, as well as in interplanetary space. In order to under- 
stand the polar atmospheric phenomena, therefore, we must know about the mag- 
netic fields around the earth and the sun. After the successful launching of arti- 
ficial satellites and subsequent discoveries of radiation belts, our old, pre-IGY 
concepts of phenomena above the ionosphere have been completely changed. 
Furthermore, many new facts have been discovered and investigated, for instance: 
the deformed earth's magnetic field holding plasma (radiation belts in the mag- 
netosphere), the outer boundaries of the magnetosphere (magnetopause), shock 
waves in front of the day-side magnetopause, and the geomagnetic tail in the 
night-side magnetopause, etc. It is impossible to describe these new environ- 
mental findings (perhaps the most exciting fields of discovery today) in this 
short review; they will be discussed at some future time. 

To indicate the present concept of the earth's environment, Figure 22 shows 
one of the latest models of the earth's magnetosphere and its tail, which has been 
built by a series of satellite observations on the magnetic field and plasma (Ness, 
1966). These models can explain complicated auroral phenomena, such as spatial 
and temporal variations of auroral arcs . However, the origin of auroral electrons and 
their acceleration mechanisms are still a great mystery. Dr. Hess of the Goddard 
Space Flight Center* (Hess, 1966) plans to produce an artificial aurora by eject- 
ing high-energy electron beams from an accelerator (order of 10 kev, pulse output 
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up to 15 kw) mounted on a satellite. This will be tested in the near future, and 
will be a great contribution to geophysics and space science. 

Finally, the importance of rocket soundings should be emphasized. Since 
such important events and phenomena as aurorae, PCA-airglow and noctilucent 
clouds, etc., are far beyond a balloon’s ceiling height (~40 km) but below satel- 
lite orbits (>100 km), rockets are the only instruments by which these phenomena 
can be observed directly. (Low-orbit, short-life satellites could be developed but 
would be less efficient and reliable.) 
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APPENDIX A 


The Ionospheric Critical Frequency 
(Plasma Frequency) 


The wave equation for the electromagnetic wave is derived from Maxwell's 
equations for Faraday's and Ampere's laws by eliminating magnetic field H. In 
vacuum this is 




c 2 Bt 2 


(A.1) 


where E is the electric field intensity, c is the speed of light in vacuum. 

In plasma, because of the mobility of ions and electrons which cause an 
electric current j* in the plasma, the wave equation is (instead of Equation A.l): 


1 3 2 E 4-rrZ pe 2 
V 2 E + — E 

c 2 3t 2 m.m c 2 

i e 


(Spitzer, 1956). 


(A.2) 


where p = n.m. + n e m e , n. and m. are number density and mass, respectively 
of ions, n and m e are number density and mass, respectively, of electrons, Z 
is charge of ion. 

Substitution of E = E 0 e i(k • r ~" t) in Equation (A.2) results in the following 
dispersion relation: 
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In the second expression of Equation A. 4, the neutrality of plasma is assumed 
(i.e., Zn. = n e ). From Equation A. 3 the phase velocity of the wave,V , in the plasma 
is given by 


(A.5) 


The index of refraction, n , of the plasma is therefore 


c 

n = — 

V 

By Equation A. 6, the index of refraction is imaginary if u> < a> . In other 
words, a wave whose frequency is lower than the local plasma frequency cannot 
propagate the medium. Since m e < < m. , the second term in Equation A. 4 can be 
neglected, and, noting f = Equation A.4 gives 

i 

f P=(“ -Y-9 • 10 3 /^T. (A.7) 

p \ m e n) 

This leads to Equation 2.1 for the F2-peak critical frequency. Generally f p , 
given by Equation A.7 is called "plasma frequency"* and is a function of electron 
concentration only. 


f- 


(A.6) 
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should be called plasma angular frequency. 
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APPENDIX B 
The Chapman Layer 


The variation of electron density, N, at a fixed altitude, z , is given by 


(IN 

dT 


q - L, 


(B.l) 


where q and L are the rates of electron production and loss, respectively. 

Loss L in Equation B.l consists of two processes: recombination of elec- 
trons with positive ions and attachment to neutral molecules. Both losses have 
rates proportional to N 2 , i.e., 


L = iN 2 , 


(B2) 


where a is a constant. 

The production rate, q , due to photo ionization of air molecules by the solar 
ultraviolet radiation, is proportional to the radiation intensity I (z), which can 
be given as a solution of the attenuation equation (Chapman, 1931), 

dl 

Y “ A ■ (z) - sec X ’ dz » (B.3) 


where A = atmospheric absorption coefficient of radiation. 

P (z) = air density at altitude z, and 
X = sun zenith angle of the sun. 

Assuming that ■ (z) / 0 -e _z/H , where H is the scale height and p Q = p(z = 0), 

the solution of Equation B.3 is 

I(z) I, i'xp -A a / q H • (sec'X) • e -z/H ), (B.4) 
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where I a is the intensity of solar radiation outside the earth's atmosphere. If 
the number of ions produced by the absorption of unit quantity of the radiation is 
y, the rate of production of ions per unit volume is 


q ( z) = y cos X — 
dz 


= AI oo/° 0 exp | -1 - Ap 0 H(secx)e z/H j 


(B.5) 


Denoting the height of the maximum production rate, by z m , which is obtained 
by dq/dz = 0, the rate of maximum ion production q m = q(z m ) is given by 


q 


m 


y 1 ^ 

H exp (1) 


cos X 


(B.6) 


where 

exp (1) y 2.71828 ... * 

The maximum number density of electrons at equilibrium is given by Equa- 
tion B.l, with dN/dt = 0 and by Equation B.2, i.e. , q m = aN 2 . 


*That is, the base of natural logarithms; its usual symbol, “e,” cannot be used here because e 
is in use as the electronic charge. 
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Combining this with Equation B.6 gives 


I 


N = 



K 

H exp(l) 



1 

cos 2 X 


(B.7)* 


Substituting this in Equation A.7 gives (note N = n e ) 

f « (cosX)^- (B.8)** 

p 

Assuming equilibrium at all altitudes gives a vertical distribution of electron 
density, which is proportional to [q(z)] 1/2 , where q (z) is given by Equation 
B.5. This is called an "equilibrium Chapman layer." 


★ ★ 


If electron loss rate L is /3N, instead of being proportional to N 2 (as in Equation B.2), the fol- 
lowing relation takes the place of Equation B.7: 


N = 



loo 

H exp (1) 


cos X, 


(B.7') 


Instead of Equation B.8: 


f « /cos X. 

p 


(B.8’) 


These equations constitute another type of “Chapman layer/’ useful in studying photo-chemical 
reactions in the upper atmosphere. The two types are called “Chapman a- layer and Chapman 
/3-layer” respectively (Ratcliffe, 1960). 
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APPENDIX C 
Ionospheric Absorptions 


The forced movements of an electron in the upper atmosphere due to the 
oscillating field of an electromagnetic wave can be described by the following 
equation of motion: 


mr = eE + er x B Q -m^r, (C.l) 

where?, E , and § 0 are vectors indicating displacement of an electron, the 
electric field of the wave, and the static magnetic field of the earth, respectively; 
m e and e are the mass and charge of an electron, respectively; and v is the average 
collision frequency between an electron and heavy particles in the medium (the 
ionosphere), called "effective collision frequency." The last term in Equation 
C.l causes electromagnetic-wave attenuation by damping the motion of electrons 
in the medium. Displacing an electron is equivalent to introducing an electric 
dipole of moment er into the medium. The electric dipole polarization, P , of the 
medium is therefore: 


P = Ner , 


(C.2) 


where N is the number density of electrons in the medium. 

Since (co is the angular frequency of the propagating electromagnetic 

wave) multiplying Equation C.l by Ne gives 

- a) 2 m P = Ne 2 E + icoeP x B - ium uP (C-3) 

e 0 c 

Following the standard notations for magnetroionic theory (Ratcliffe, 1959), this 
equation becomes 



XE = P(l-iZ)-iY x P, 


(C .4) * 


*ln Gaussian units, the electric permittivity of free space e Q is unity. 
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where 


X = 


477Ne 2 


a > 1 


e 0 m e ar \ co 


(cf. Equation A.4) 


Y = - 


e„B, 
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and 


Z :=Ji 


Taking the Cartesian coordinates x ,y , z (where z\ |k , k is propagation 
vector or direction of wave normal), y| | B 0 , and writing Y L =Y cosd , Y T = Y sin 6, 
6 is the angle between k and B 0 ) gives Equation C.4 as 
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(C.5) 


Since all components of the wave fields vary as e ikz (a plane wave prop- 
agating along the z-axis), Maxwell's third and fourth equations 


,0 ® ? 3H 

curl H ; - — , curl E ; - u 

0 Bt 


Bt ’ 


can be written as follows (noting thatD = e 0 E + 4 -ttP , where m 0 ; magnetic 
permittivity of free space): 


- ikH = Leo (e n E + 4?rP ) , 

x v Q v v 


ikH y = Leo (e 0 E x + 4 ttP x ) , 
0 = Leo (e 0 E z + 4 ttP z ) ’ 


(C.6) 
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and 


- i"M 0 H x - ikE y - 


i^ 0 H y = ikE x , 


i"M 0 H z = 0 • 


Equations C.6 and C.7 give 
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Since the phase velocity of the wave isV =oj/k = c/n , and because < 
Equation C.8 becomes 

2 - P x 

n 2 = 1 -r~. 

e o ® 

0 x 


From Equation C.5 


o r 


can be eliminated in Equation C.9, giving 


n 2 = 1 - 


1 - iZ - iY L R 


where 



Since the index of refraction n is complex, this can be written 


(C.7) 

(C.8) 

~ e o * 
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n = /j. - it] • 


(C.12) 



Then, the wave attenuation is e*^ after propagating a distance z = c/u = k/2i 
with velocity c ■' ji. In other words, the absorption co-efficient k, defined by 


E z; E 0 e“/<z 


C 0 

K • T). 

c 


is given bj' 

Combining Equations C.10 and C.12 gives 

1 47iNe 2 
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Another expression is 
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At places where fx»l, Equation C.16 can be written, assuming 2 < < 1 ~. 2 , 


2c (1±|YJ) 


A 


(C.17) 


In the above expressions, the plus sign represents the ordinary wave; the 
minus sign represents the extraordinary wave. 

In most of the ionosphere, it can be assumed that v 2 << ( a. + UJ ) 2 and, if 
o;>>a; L , (i.e., Y<<1). p- - 1 then 
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The reflection coefficient, p , can be given approximately as 
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where N = N (z), v = v(z), and hj is the height of total reflection of the wave 
with angular frequency 03 . 


An absorption obeying this frequency law (i.e., Equations C. 18 and C.19) 
is called a "non-deviative" absorption. 

On the other hand, when Y L < 1, Equation C.17 can be written: 



which gives, instead of Equation C.19, 


where 
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(C.21) 


are called equivalent (or group) height, and optical (or phase) height, respectively. 
Absorptions expressed by these terms are called "deviative" absorptions. It 
should be noticed that in the case of non-deviative absorption /^-l ; therefore 
h - h r and/< ~ 0. Deviative absorptions are negligible; only non-deviative 

K * 

absorption is important in frequency ranges -1. 
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APPENDIX D 

The Units of Auroral Brightness 


Conversion of wave length A (in A) of the light emitted in the transition be- 
tween two energy states E. , E f (in eV ) of atoms or molecules is obtained from 
the Planck's photo quantum formula 


hv = AE, 


(D.l) 


where h = 6.625 x 10 27 erg see(~4.15 x 10 -15 eV sec), 


£77 


= — • c (c is light velocity, c ^ 3 x 10 10 cm/sec) 

A. 


and AE = E . - E f . 


Equation D.l is therefore written 


A.AE = 1.24 x 10 4 (eV A) 


(D.2) 


For example, the energy of primary electrons AE required to excite oxygen 
atoms to the state that emits red light, 6300 A (see Type A aurora) is, from 
Equation D.2, 


AE - 2 eV • 


This is indicated in the energy diagram in Figure 11(a). 

The brightness of aurorae and general airglows are generally expressed by 
the units called Rayleigh (R), defined by 10 6 photons per cm per sec. This unit 
is based on the total emitted photons in the column along the line of sight toward 
the source per unit solid angle (steradian) at the observer, multiplied by 477 
(i.e., assumed isotropic luminosity 
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where F(r) is rate of emission at distance r from the observer (photon/cm 3 
sec); the value of 477 J is "one Rayleigh" if it is 10 6 in photons per cm2 per sec. 
In other words (Chamberlain, 1961), 1 R = an apparent emission rate of one 
megaphoton per cm2 per second. 

As can be seen from Planck's formula (D.l), conversion to energy flux from 
Rayleigh units differs for different wave lengths of emitted light. Unless other- 
wise specified, the brightness of aurora is based on the intensity of oxygen line 
5577 A, which has been used to define the International Brightness Coefficient 
(IBC): IBC I, II, HI, and IV corresponds to brightness (10 3 R, 10 4 R, 10 s R, and 
10 6 R respectively. 

The energy flux of light emitted by the IBC-I (i.e., weak) aurora is, by the 
quantum Equation D.l, 3.6 x 10 - 3 - 4 x 10" 3 ergs/cm 2 sec. 


Table D.l 

Energy Fluxes Around the Earth 


Source 

Energy flux (ergs cm -2 sec -1 ) 

Sun 

1.4 x 10 6 (solar constant) 

Full moon 

3 

Total starlight 

1.8 x 10- 3 

Airglow (visible) 

1.5 x 10" 2 

OH-airglow (infrared) 

2 x 10" 2 

Lyman alpha (ultraviolet) 

10-2 ** 

Cosmic rays 

3.8 x 10~ 3 


*Comparable with the brightest (IBC-IV) aurora. 

**The solar Lyman Alpha emission (approximately 5 ergs cm'^ sec -1 ) is excluded; because of 
strong absorption by air, this is not observable from the ground. 
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Figure 1-Diurnal variations of f Q F2 in the Antarctica based on the monthly medians of hourly 

values (December, 1957). (Hill, 1960). 
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igure 3-Effects of wiTid shear and geomag 
netic field in the F-region. 



cosx 


Figure 4-Average monthly values of absorption versus inverse Chapman’s function 
[ch cos for constant sunspot numbered R - 100 and 200, respectively 

(f -- 4 MHz). 
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Figure 7-Local time-dependence of electron precipitation (E > 40 kev) by satellites 
Injun 3 (Frank/ Van Allen, Craven, 1964). 
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Figure 8-Temporal and spatial (horizontal) distribution of the 
auroral display observed at College, Alaska, on January 
25-26, 1952 (Heppner, 1958). 
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RELATIVE INTENSITY, KR ABSORPTION, DB 



Figure 13-Time variation of polar cap glow aurora during 7-23, July 1959 reported by Sanford 
(1962), compared with K-index of Scott base in Antarctica, Solar flare and planetary K p index, 
Thule-riometer records (Am 88°N) and intensities of zenith auroral emission (in KR); thin line 
5566& [01] , and thick line 3914& [N 2 + l line, observed at Scott station in Antarctica (79 C S). 
It should be noted that riometer data is from the northern hemisphere, aurora data from the 
southern. (Sanford, 1962). 
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* MT MCKINLEY 



19 FEB 1966 


Figure 15— Horizontal projection of a fast-moving auroral arc at one-minute intervals, observed 
at Col I ege, Alaska, on February 19, 1966. (Time attached to arc is AST - UT - 10). (Wilson and 
Nichparenko, 1967). 
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Figure 17-Observed and calculated polarization (percentage polarization, P) of blue light from 
noctilucent cloud versus scattering angle i//. Full line corresponds to the Rayleigh scattering; 
dashed lines are calculated results based on Mie*s theory. Four values of parameter x: x “1.2 
and 1.5 correspond to particle-radius r = 0.09 /x and r = 0. 1 1 jJL respectively, with n = 1.33 (water); 
x -(1.5) and (2.0) correspond to r~0.ll/LX and 0.16/x with n = 1.55, respectively, n is the 
refractive index of particle. Two thin lines at the bottom indicate the magnitude of 5% error.(Witt, 





Figure 18-Geometric conditions illustrating conditions for 
seeing the noctilucent cloud (NLC) from the ground, where 
a is the solar depression angle, and H sc is atmospheric 
screening height. (Fogle, 1966). 
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JANUARY FEBRUARY MARCH APRIL MAY JUNE JULY AUGUST SEPTEMBER OCTOBER NOVEMBER DECEMBER 
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NOCTILUCENT CLOUD OBSERVATION PERIODS {SHADED AREA) FOR 40*N. 


Figure 19(a)-The observable period of noctilucent cloud at different 
(geographic) latitudes, 40°. 
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Figure 19(b) — The observable period of noctilucent cloud at different 
(geographic) latitudes, 60°. 
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Figure 20-Vertical temperature 
distributions observed by rocket 
(grenade) soundings at Ft. 
Churchill (A.-59°N, 68 ° 6N) 

Canada. (Stroud et al., 1960). 
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Figure 21 —The stratospheric diurnal circulation at the sum- 
mer sol stice, projected on the equatorial plane. Thin solid 
curves are contours of constant temperature or thickness; 
the arrows indicate the tidal wind field. This diurnal circu- 
lation is imposed on the general summer easterly and weak 

r i rn i I nt! An fWakk 10AA 






